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We construct and probe a holographic description of state of matter which results from coupling 
a Fermi liquid to a relativistic conformal field theory (CFT). The bulk solution is described by 
a quantum gas of fermions supported from collapse into the gravitational well of AdS by their 
own electrostatic repulsion. In the probe limit studied here, the Landau quasiparticles survive this 
coupling to a CFT. 



INTRODUCTION 

Some progress has been made recently in using holog- 
raphy to learn about metallic states of matter which are 
not described by Landau's Fermi liquid theory [1-14]. 
In particular, the work of [2, 4] constructs controlled 
non-Fermi liquid fixed points, which however constitute 
a parametrically-small (in the bulk Newton's constant 
l/N"^) fraction of a larger system; this larger system is 
responsible for the destruction of the Landau quasiparti- 
cle. 

The fact that this larger system is locally critical is the 
origin of the short transport lifetime and linear-T resistiv- 
ity in the case which realizes marginal Fermi liquid [5, 6]. 
This local criticality is further closely tied to the nonzero 
zero-temperature entropy of this state, indicating that in- 
stabilities prevent one from cooling into this state. Var- 
ious instabilities have been suggested; an intrinsic one, 
arising from the density of fermions itself, was pointed 
out in [8] : these fermions screen the gauge flux which sup- 
ports the ^^5*2 throat, and a truer groundstate is a Lif- 
shitz geometry with large dynamical exponent z ^ 7V^, 
rather than the 2; = oo of AdS2 x M^. Such 'electron 
star' states, comprising a charged and gravitating den- 
sity of fermions in the bulk, have been further studied 
in detail [9-14]. These states represent an improvement 
over the work of [2, 4] in that the fermions play a leading 
order (in N"^) role in constructing the geometry. How- 
ever, the single-fermion response in these states exhibits 
many Fermi surfaces [12, 14], rather than one. This con- 
clusion is not intrinsic to holographic states supported by 
fermions, but rather is an artifact of the Thomas-Fermi 
approximation used in the pioneering work on electron 
stars [8-14]. 

We would like to construct a holographic non-Fermi 
liquid (NFL) where the fermions contribute at leading or- 
der to the state {e.g. to the thermodynamics and ground- 
state entanglement entropy of regions). 

A point of departure is provided by [15]^. The state 
constructed there can be understood as a Fermi liquid 



^ Earlier work which studies quantum spinor fields in a holographic 
context includes [16, 17]. Recent provocative work towards this 
goal includes [18—21]. 



coupled to a (toy model of a) confining gauge theory. The 
outcome is a Fermi liquid, in bulk and boundary. The 
lifetimes of the quasiparticles are infinite in the leading 
large- limit. 

This paper makes a further step towards the above 
goal for holographic NFLs. We construct a bulk Fermi 
liquid state without putting in a gap in the geometry 
from the beginning. That is, we show that the bulk state 
constructed in [15] survives the limit where the hard-wall 
cutoff is taken away. The problem we solve is the fermion 
analog of the probe holographic superconductor calcula- 
tion in [22]. The resulting state is a quantum electron 
sta'T in the sense that we are treating the bulk fermions 
quantum mechanically, but they do not collapse into the 
Poincare horizon. Solving this problem is a prerequisite 
for the more general problem including gravitation. 

Technically, this required some improvements over pre- 
vious methods. In particular, in order to understand ac- 
curately the effect of the bulk electric field on the charge 
of the filled Dirac sea, it was necessary to provide a short- 
distance (UV) completion of the bulk system. We accom- 
plish this by putting the bulk system on a lattice. As a 
result, we were forced to account for the polarization of 
the Dirac vacuum by the electrostatic potential, and the 
resulting screening of the electric charge. 

Next, we describe the problem, and outline our method 
of solution, deferring details to the appendices. Results 
and discussion follow. The appendices also include a dis- 
cussion of several interesting effects uncovered here, in- 
cluding the chiral anomaly in the bulk, and a surface 
charge density. 



SETUP OF THE PROBLEM 



We consider a system defined by the action 



S = / d-'+'xy^ 

where 



R-2A 1 



IGttGn 4g2 



= / e 



iSf 14>,A] 



+ W[A] , (1) 
(2) 



This usage has appeared in [12]. 
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tp = V'^T^ and Vm = du + j'^abA/i 
the spin connection. 

We study this system in a probe hmit, Gn — > at 
fixed A, where the geometry is not dynamical. In the 
dual language, we are studying a CFT where the fraction 
of degrees of freedom that carry charge is small. 

In particular, we specialize to the AdS4 metric 

ds^ = -L (_di2 + dx^ + dy2 + dz2) , (4) 

and we consider a gauge field of the form A — ^{z)dt. 
The equation of motion for the potential <& is 



(5) 



To compute the expectation value on the right hand 
side, we we expand the spinor field operator in eigen- 
functions of the Dirac Hamiltonian. The Dirac equation 
is {If) + m)il) = 0, and we make the following ansatz for 
the eigenf unction tp: 



^ = z3/^e-'"^^+''^-'^^'{z) . 



(6) 



For any fixed fc, the Dirac equation can be block- 
diagonalized [4]. Without loss of generality we set k — 
kx, and we let ^1/ = (5"+, ^_)*. Using the following basis 
for the Clifford algebra: 



we have 
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Eq. (8) has two linearly independent solutions, whose 
asymptotic behavior near the AdS boundary is 



(9) 



plus terms subleading in z. We demand that the non- 
normalizable solution be zero^: a = 0. Moreover, as in 
[15], we impose hard- wall boundary conditions at an IR 
cutoff z = Zm'. the upper component of both '^±(zm) 
must vanish. 

With these boundary conditions, the Dirac Hamilto- 
nian is a self-adjoint differential operator with spectrum 
'^n fe s' l^^belled by a discrete index n, by the momentum 
k and by the sign s = ± that distinguishes the upper 



components from the lower components. We denote the 
eigenfunctions with 5* r ■ 

o n,k.s 

In order to give definite meaning to the expectation 
value in (5), for a given profile of of <&, we fill all the states 
with LOn_k,s < 0, and we subtract the same expectation 
value for $ = 0. That is, we solve 



$"(z) = q2 [^(^)i^ _ n{z)\^=o\ = q^An{z) , (10) 



with 



n{z) = J2 J*l,^,,(^)*a^.(^) • (11) 

Clearly (11) involves two sums that need to be regu- 
lated. We regulate the sum over n by discretizing the 
z coordinate, with lattice spacing Az, and we impose a 
hard cutoff fc^ < A| on the sum over momenta. After 
the subtraction in (10), and appropriate renormalization 
of the charge q, the problem has a well-defined limit as 
Az — > 0, Afc — > oo. Additional information on the sub- 
traction and renormalization can be found in the supple- 
mentary material. 

Eq (10) also needs to be complemented with appro- 
priate boundary conditions on $. We want a finite 
chemical potential in the boundary theory, so we set 
$(0) = —fi, and we also impose ^'(zm) — 0. We expect 
that the boundary condition at z„i becomes unimpor- 
tant as Zm and we verified this by exploring mixed 
boundary conditions as well, without finding a significant 
influence in the interior, for large enough Zm- 

We solve the integro-differential system formed by (8) 
and (10) by an iterative method, whereby the number 
density computed with a given a profile of $ is used to 
update $ through (10) and the new profile is used to 
update the number density, until convergence. 



THE RESULTING GROUNDSTATE 

Fig. 1 and 2 display typical profiles for $(z) and n{z). 
They also show that the problem possesses a well-defined 
limit as z,„ — )> cxD, Afc — > oo. We also verified that the 
profiles are insensitive to the discretization of the z coor- 
dinate, for small enough lattice spacing. Once the cutoffs 
are removed, the only remaining scale in the problem is 
the chemical potential /x, and, without loss of generality, 
we can set /j, = 1. 

The profile of the potential approaches a constant as 
z — oo"', and we would like to argue that this constant 
is zero. 



^ Allowing the mass to range over (— |,oo), this includes the al- 
ternative quantization. 



* If the electric field oc <&' does not vanish at the Poincare horizon, 
an argument similar to §7.4 of [8] indicates that there will be 
a value of z beyond which backreaction cannot be ignored, no 
matter how small the Newton's constant. 
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FIG. 1. This plot demonstrates the existence of the hmit 
— >■ oo with, from red to blue, Zm. ~ 10, 20, 30, 40, 50. 
The associated values of the boundary chage density are 
p = 0.1710, 0.1669, 0.1638, 0.1627, 0.1622. m = 0.3; q = 2.0; 
dz = 0.1; Afe = 20. 




FIG. 2. This plot displays how the limit Afe — >■ oo is ap- 
proached, with, from red to blue, — 2, 3, 5, 10, 20. See 
the appendix for further discussion of the surface charge den- 
sity. Corresponding values of the boundary charge density 
are p = 0.1721, 0.1717, 0.1714, 0.1710, 0.1708 (note that here 
Zm = 20, so these numbers should be compared with the 
first value in fig. 1). Positive k represent the spectrum of the 
m = 0.3; q = 2.0; Zm = 20; dz = 0.04. 

The asymptotic behavior of the wavefunction for large 
z depends on the sign of 

q2 = (c^ - $(oo))2 ~ fc2 _ (^2) 

For < the wavefunction is exponentially decaying, 
whereas for > 0, ^'(z) ^ e*''^. Let us call the region 
q^ > the infrared (IR) lightcone. As we increase Zm, the 
gap between the bands inside IR light cone decreases like 
l/z„i. Hence, for — >■ oo, a continuum develops inside 
the IR light cone (fig. 3). The contribution to the number 
density coming from each state within the IR light cone 




FIG. 3. The portion of the spectrum of the Dirac operator 
nearest to the chemical potential. For fc > we display LJn,k.+ , 
for fc < we display tOn^k,-- Also shown are the IR lighcone 
(dashed lines) and the UV lightcone (dash-dotted lines). It 
is clear from the figure that a continuum is developing inside 
the IR lightcone. m = 0.3; q = 2.0; Zm = 200; dz = 0.2; 
Afe = 20. 



also decreases like 1/zm, and is constant in z, for large 
z. Therefore, any finite portion of the lightcone that lies 
below uj = gives a finite, z-independcnt contribution to 
the number density. 

If $(oo) < 0, after the subtraction of (10), the lower 
half of the light cone does not contribute any net num- 
ber density, but the portion of the upper half of the light 
cone that lies below uj — gives, at large z, a finite, z- 
independent contribution to the number density. This is 
incompatible with the potential going to a constant as 
z — >■ oo, due to Gauss law. A similar argument obtains 
for the case $(oo) > 0. The only possibility is for the 
potential to go to zero, so that the states within the IR 
light cone do not contribute at all to the number den- 
sity density. Then, the only contribution comes from the 
few bands that lie outside the IR light cone. The corre- 
sponding wavefunctions are exponentially decaying, they 
don't contribute to the number density at large z, and 
hence they are compatible with the potential approach- 
ing a constant (zero) as z — )> cxi. 

Determining numerically the exact nature of the falloff 
of the number density at large z is very difficult, but it is 
quite manifestly subexponential. A possible explanation 
for this is that the bands that lie outside of the IR light 
cone, for a certain range of k, skirt the edge of the cone 
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(follow the dashed line), as can be seen in fig. 3. When 
the state is close to the edge of the cone, the rate of 
decay is very weak. Since the distance from the edge is a 
decreasing function of Zm, it is conceivable to obtain a a 
subexponential decay of the number density as — >■ oo. 

From the holographic point of view, the IR region of 
the geometry is dual to a relativistic CFT, which only 
has spectral weight inside the lightcone |w| < ck. One 
interesting quantity that can be extracted from our com- 
putation is the boundary charge density p. This is the 
response to the chemical potential, and it is given by 

p = <!>'{0)^q^ J dz An{z), (13) 

where the second equality is a consequence of Gauss law. 
The charge density in the boundary is equal to the total 
charge in the bulk. Since p is the only scale in the prob- 
lem, its dependence on /i is determined by dimensional 
analysis to be p = Ap,^, for some constant A{q, m). 

Let us also point out that (13) guarantees Luttinger's 
theorem in the boundary [15]. Luttinger's theorem states 
that, for interacting fermions, the area of the Fermi sur- 
face is proportional to the number density, with the same 
factor as for free fermions. As will be discussed in more 
detail in the next section, there is a Fermi surface wher- 
ever one of the bands in fig. 3 crosses a; = 0. According 
to our construction of An(z), each k mode within the 
Fermi surface contributes 1 to integral on the RHS of 
(13), thereby ensuring Luttinger's theorem. 

GREEN'S FUNCTIONS 

We compute the photoemission response of our 
state, proportional to the single-fermion spectral den- 
sity lmG}i{uj, k), where Gh is the retarded single-fermion 
Green's function. To compute the retarded function, we 
impose in-falling boundary condition for the Dirac field 
at the Poincare horizon, and include a source at the UV 
boundary in (9), a ^ 0. Then Gr{uj, k) = h/a. See fig. 4. 

The lifetimes of quasiparticles in holographic Fermi 
surfaces can be understood in terms of interactions with 
other gapless degrees of freedom [4, 6, 7, 23]. The de- 
grees of freedom into which a Fermi surface quasiparticle 
might decay are those inside the IR lightcone described 
above. 

Outside of the IR light cone, for uj = ujn,k, there ex- 
ists a solution to the Dirac equation that is real, decays 
exponentially at large z and is normalizable in the UV. 
This means that the infalling boundary condition is triv- 
ially satisfied, because the wavefunction is zero at the 
horizon, and we have a finite response 5 7^ for zero a. 
As a consequence, Gu^uj, k) has a pole dX uj — uJn,k, and 
hence there is a delta function singularity in ImG^. This 
delta function implies the existence of an exactly stable 
quasiparticle in the boundary theory. In particular, the 
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FIG. 4. Density plot of the spectral density ImGij(cj, fc), 
displaying the spectra of the stable quasiparticles, which co- 
incide with the bands in fig. 3, and the continuum inside the 
IR light cone. Notice the increased width of the quasiparti- 
cle peaks as they enter the light cone, m = 0.3; q = 2.0; 
= 200; dz = 0.2; = 20. 

points where — are Fermi surfaces with stable ex- 
citations. To detect such infinitely-narrow resonances in 
the numerics, we add a small imaginary part to the fre- 
quency, so as to move the pole away from the Re(aj) axis, 
and convert the delta function to a narrow Lorentzian. 

On the other hand, inside the IR lightcone, the asymp- 
totic behavior of the wavefunction is 4* ^ e*^', and 
hence, in general, a finite and complex Gij(a;, fc) is needed 
to satisfy infalling boundary conditions. Consequently, 
quasiparticle excitations have finite width, because they 
can decay into the gapless CFT excitations. This phe- 
nomenon is visible in fig. 4 where the bound state bands 
enter the IR lightcone. 



DISCUSSION 

The state we have discussed can be described semi- 
holographically [4, 7] as arising from a Fermi liquid cou- 
pled to a relativistic CFT. The study of Fermi surfaces 
coupled to critical systems has a long history, e.g. [24- 
31]. The conclusion of our holographic calculation is that 
the coupling between these sectors described here is an 
irrelevant deformation of the Landau theory. In fact, ac- 
cording to our discussion about the stability of the quasi- 
particles, and the location of the IR light cone, the only 
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possible singularities at w = 0, 7^ are delta func- 
tion peaks, which indicate exactly stable quasiparticles. 
This fact is very likely a consequence of the probe limit 
Gn 0. 

The nature of the coupling between the FS and the 
CFT that one infers for the semi-holographic picture 
is a hybridization between a fermionic operator of the 
CFT and the electron operator, as in [4, 7]. Possibly- 
relevant couplings between the fermion density and rele- 
vant bosonic operators of the CFT of the kind considered 
in [32] are suppressed in our large- limit. 

It will be very interesting to study the effect of the 
screening by the fermions on the boundary gauge the- 
ory dynamics. On general grounds we expect that, with 
enough fermions, even beginning with a confining solu- 
tion a,t n — 0, the gauge theory will deconfine. Holo- 
graphically, this requires taking into account the gravi- 
tational backreaction of the bulk fermions. Progress in 
this direction will be reported elsewhere. Resolving the 
problem confronted in this paper - the question of the 
state of the bulk fermions in the presence of a horizon - 
was an essential prior step. 
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SUPPLEMENTARY MATERIAL 



Renormalization prescription 

Here we describe our method to obtain a number den- 
sity which is independent of the cutoffs Ak and Az, and 
whose z integral is also cutoff-independent. 

Since the ultraviolet divergences probe short-distance 
structure of the theory, we can understand them while 
working in the flat-space limit, i.e. to lowest order in an 
expansion in curvature. The expectation value of the 
fermion number current in a given background can 
be evaluated in perturbation theory as 
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(■07''-0) = ^^^^ + ^^^^^^^^)ww + {^^^^l)"'^ + ■■• ' 

where the box represents an insertion of the current, and 
the wiggly hnes represent the background field. Count- 
ing powers of loop momenta, diagrams with up to four 
insertions may be divergent. The first term is made 
finite by the subtraction of the zero-field value as in 
(10); this ensures that the density in vacuum is zero: 
Ans = ~ (■!/'^'0)vi=o- Following familiar analy- 

sis from QED in flat space, loops with odd numbers of 
insertions vanish in the UV limit by charge conjugation 
invariance, and the divergence in a loop with four inser- 
tions is ruled out by the Ward identity. In the remaining 
loop with two insertions, the naively expected quadratic 
divergence is again ruled out by the Ward identity, and 
the only remaining divergence is logarithmic. Specifi- 
cally, the bare density for a background field of the form 
A = ^{z)dt is 



AnB(z) = 



5W[A = ^dt] 
(5$(z) 
d-^n ~ 

— $(p)e-*P^^; 



(2^)4 



tr [7"L)F(fc)70^f (A: +p)] + 0($2) 



dHE 1 
{2^Y 1% 



^"(z) -I- finite terms 



(14) 



where ^{p) has delta-function support at Pt — Px ^ Py = 
0, Dp is the momentum-space fermion propagator, and 
the subscript B denotes a bare quantity. 

Thus, with a Lorentz-invariant cutoff A, the logarith- 
mic divergence in the bare density is 



Ansiz) 



1 

(3^ 



In ^"{z) -{- finite terms , (15) 



where M denotes some infrared energy scale. In our non- 
relativistic numerical setup with two cutoff scales A^, 
1/Az, the divergence has a more complicated form 



AnB{z) = x(Afc, Az) <l>"(z) + finite terms , 



(16) 



where x diverges logarithmically with s under the simul- 
taneous scaling A^ — > sA^, Az — > Az/s. In the next sub- 
section we use perturbation theory to examine in more 
detail the number density and the function x, in the limit 
Az 0, Afe finite. 

In the numerical analysis, we compute x(Afc,Az) 
at a given value of the cutoffs from the charge re- 
sponse to a small background field of conventional pro- 
file ^i,{z) = ez^ . In order to be able to neglect the fi- 
nite, nonlinear terms in the response (16), we require 



£ <C z'^\ the quadratic profile conveniently removes any 
higher-derivative contributions. 

Once X is known, for any given we define a renor- 
malized number density as 

Anj,{z) = AnB{z) - x(Afc, Az) [$"(z) + <5(z)<i>'(z)] . 

In the iteration procedure, the electrostatic potential for 
the next step $ncxt is the potential produced by the 
renormalized density according to 



<cxt(2) -4,Anfl(2) 



(18) 



The renormalized charge q^j is by definition cutoff- 
independent - it is the charge parameter specified in the 
plots above. Apart from the delta function term, the sub- 
traction (17) has the effect of renormalizing the charge q 
according to 



4 = 



l + 'z|x(Afe,Az) ' 



(19) 



as can be readily seen by substituting this and (17) into 
(18). Away from the boundary, this implies q^Anji — 
q^AuB- 

The delta- function term in (17) is required by the fol- 
lowing logic. Its purpose is to guarantee that integrating 
both sides of (17) over z gives the relation 



J dz Anji{z) = y dz AnB{z) 



(20) 



This is very important, because the right hand side does 
not depend on the cutoff, for large enough A^ (larger 
than any Fermi momenta, so that varying A^ does not 
change which states are filled relative to the vacuum), 
since the contribution of each additional fc-mode to the 
integral is then exactly zero. On the left hand side, the 
renormalized charge density should also not depend on 
the cutoff, hence the need to add the delta-function term 
at the boundary. 



Single mode linear response 

In this section we compute at first order in perturba- 
tion theory the number density response of a single fc- 
mode, in the limit Az — > 0. This is important as a check 
on the numerics, and for the interpretation of the surface 
charge in the next section. We will assume translational 
invariance, so our results will be valid when fc 3> 1/zm, 
for ^ z <C Zm- In this regime the bulk fermion mass 
has no effect, so we focus on m = 0. 

We define the susceptibility x as 



Anfe(z) = / d7Jx{z,z'Mz'), 



(21) 
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and, at first order, we have: 



m|wT,^>0 n|cj„<0 



where labels a basis of eigenstates of the unperturbed 
single particle Hamiltonian h — ka^ — ia^dz- The eigen- 
modes are 



with LOq E 

tors are 



■^qs , {kcr^ + gcr^)*gs = SW^^gs , (23) 



^y q'^ + k'^ and s = ±1. The explicit eigenvec 
1 



y^2ujg{LUq ~ k) ~ ^ 

1 



UJq — k 



yj2ujq{iUq - fc) V ^ 

Plugging these expressions into (22) gives 
X(^,^') = -|^ / ^e^'^(^-^'^T(p)+h.c. 



(24) 
(25) 

(26) 



with 



4arctanh 



Tip) = 1 



(27) 



The fourier transform of T does not exist, but we can 
introduce a different susceptibility x 



Anfc(z)= / Az'x{z-z'W{z'). 



h.c. 



and we have a well defined answer 



X{z) 



27r 



27rg2 



(28) 



(29) 



Using the explicit expression (27) for T, we observe 
numerically that 



X{z) w 

ak 



-b\kz\ 



aw 18.16,6^2.17, 



(30) 



that is, the response of a mode of momentum k is expo- 
nentially localized to a width of order in the radial 
direction. 

We verified that the numerical computation of the 
number density, at large fc, matches this result down to 
the coefficient a and h. In particular, figure 5 shows the 
number density response to a source such that <i>"(z) = 
5{z - z™/2). 

In the special case of = ^z^, which is (proportional 
to) the conventional profile we use to compute x(Afc, Az), 
we have 



dz x{z) 



1 T(p) 
— — hm — — 



1 



Qnk'^ 



(31) 



0.004 



0.002 



(22) ^ 



-0.004 




FIG. 5. In this figure we show the contribution to the charge 
density from modes with various momentum fc, in response to 
a profile of the electrostatic potential which has delta-function 
in the second derivative. The profile of "1> is indicated in the 
inset. Notice how the response to the source is similar to 
the finite size effects at the boundary, k from 10 to 100 in 
logarithmic steps. 



Integrating this result over momenta reproduces the fa- 
miliar logarithmic divergence (15) [q.v. [33], eqn (7.96)). 
More explicitly we have: 



x(Afc,Az = 0) = 2 



cPk 1 



1 , Afc 
log 



M 



(27r)2 67rfc2 67r2 " M 



(32) 

where M is again some infrared scale, and the factor of 
2 comes from the sum over spins. 



The surface charge 

The bulk system we are studying can effectively be 
compared to a metal. It possesses an (infinite) num- 
ber of completely filled bands, which we will call valence 
bands, and a number of partially filled bands, which we 
will call conduction bands. If a transverse electric field 
Ei is applied, there will be a current response. How- 
ever, the boundary conditions prevent the charge from 
moving in the z direction, like in an isolated slab of 
metal. Therefore, the fermions respond to the electric 
field Ez — — <I>'(z) < by attempting to screen it, in a 
way that can be compared to the polarization of an insu- 
lator. This is the origin of the thin layer of charge that 
is visible close to the AdS boundary in figs. 1,2,5,6. 

The contributions to the layer from the valence and 
the conduction bands is very different, as can be seen in 
fig. 6. In particular, the contribution from the valence 
bands forms a thin layer, which is largely independent of 
the fermion mass; we infer that it is sustained by wave- 
function gradient pressure, i.e. the uncertainty principle. 
On the other hand, the contribution from the conduc- 
tion bands depends strongly on the fermion mass, and a 
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FIG. 6. This plot shows the separate contributions to the 
charge density from the conduction bands (sohd line) and 
valence bands (dashed line), with, from red to blue, m = 
—0.3, 0.0, 0.3, 0.6. The corresponding values of the boundary 
charge density are p = 0.3105,0.2165,0.1623,0.1304. q = 2.0; 
= 50; dz = 0.071; Afc = 20. 



positive mass has the effect of pushing away the charge 
density from the boundary. 

From the point of view of the field theory dual of our 
bulk system, this surface charge at the boundary of AdS 
can be interpreted as the screening due to the rearrange- 
ment of the short-wavelength modes of the CFT. 




5 10 15 20 



FIG. 7. This plot shows the dependence on q, with, from red 
to blue, q — 0.5, 1.0, 1.5, 2.0. The corresponding values of the 
total charge are p = 0.0584,0.1055,0.1411,0.1623. m = 0.3; 
= 50; dz = 0.071; Ak = 20. 



Fig. 7 shows the effect of the charge on the total num- 
ber density and on the profile of the field. A larger charge 
leads to a more rapid screening of the chemical potential. 
The surface charge also becomes comparatively larger, 
due to the stronger electric field at the boundary. 
The bulk chiral anomaly 

An important check of our regularization of the 
Dirac Hamiltonian is the ability to reproduce the chiral 
anomaly, in the following way. 

Using the Dirac equation in the continuum, one can 
show that, for each mode 'i>n.k,s = ^ and for any choice 
of <&, the following equation holds: 

= 2*^ (-mcr^ ± /ccr^) * . (33) 

A naive application of this statement would tell us that, 
for m = and fc = 0, the profile of the charge density is 
a constant in z: 

d^nk = for m = fc = 0. (34) 

This statement is in fact false. Indeed, the unregulated 
statement (34) asserts the conservation, at the classical 
level, of an axial current. This conservation law is vio- 
lated quantum mechanically, in the presence of a back- 
ground gauge field coupling to the vector-like current. In 
fact, if we limit ourselves to study the k — mode, the 
bulk Dirac problem is that of a 1+1 dimensional fermion 
(or rather two such fermions, because the spin index be- 
comes a species index), described by the action 

S = J dzdt i^ip + m)^ (35) 

with p = + 7* {dt + i$). 

When TO = 0, this action is invariant under chiral 
rotations ^' i— >■ e*"'*' and hence the axial current 
= 'I'^^^'^v]/ is classically conserved. If we define 
7^ = 7^7*, we have jf = ^f^^l/^ and hence (34) states 
the conservation of the axial current. 

At the quantum level, the axial current conservation is 
violated by the chiral anomaly: 

d,ji: - ^e.^F'^- ■ (36) 
If we integrate this equation with respect to z we get 

nk{z) ^ — — ^ + const for m = k = 0, (37) 

TT 

which we observe numerically. 



